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Abstract 



We study a highly disordered network of superconducting granules linked 
by weak Josephson junctions in magnetic field and develop a mean field theory 
for this problem. The diamagnetic response to a slow variations of magnetic 
field is found to be analogous to the response of a type-II superconductor with 
extremely strong pinning. We calculate an effective penetration depth X g and 
critical current j c and find that both A" 1 and j c are non-zero but are strongly 
suppressed by frustration. 
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In the physics of type II superconductors it is common to distinguish between weakly 
and strongly disordered limits of the mixed state. The former case is described by the flux 
lattice distorted by disorder, while in the latter case the vortex loops proliferate and more 
adequate description is provided by the model in which the set of granules is coupled by 
Josephson junctions. Whereas the weak disorder limit has been discussed extensively [|TJ, 
the opposite limit received much less attention. It is not a priori clear how these limiting 
cases are related. In this paper we apply the methods developed in the theory of spin 
glasses and show that response of the Josephson network to the variations of magnetic field 
is qualitatively similar to the response in a weakly disordered limit. The goal of this paper is 
two-fold: it derives the observables characterizing the Josephson network in strong magnetic 
field (effective penetration depth A g , critical current j c and critical field H c ) and it establishes 
parallels between phenomena in a vortex states (depinning current, shielding, formation of 
Bean state) and phenomena in the Josephson networks. 

The previous treatment of frustrated Josephson arrays 0, [[| did not consider the spa- 
tial (or temporal) variations of magnetic field making impossible any comparison with 
the experiment which typically measures the magnetic moments induced by screening 
currents. Here we extend the MF theory so that it includes the magnetic field varia- 
tions, the resulting MF equations are analogs of Ginzburg-Landau equations for the glassy 
state. Unlike the Ginzburg-Landau equations which are exact in a wide temperature range 
Gi <C {T c — T)/T <C 1 (where Ginzburg parameter Gi <C 1), the MF equations for the glassy 
state hold only if the effective range £o of the Josephson couplings is large compared to the 
distance r between granules. This condition is not realized in conventional granular super- 
conductors, but we hope that results obtained in this approach are still qualitatively correct 
and provide reasonable quantitative estimates |3J]. When deriving the MF equations we shall 
also assume that the external magnetic field varies slowly compared to a typical relaxation 
time of each junction and that its spatial variations are small on the scale of £o and that 
individual loops of Josephson links are too weak to keep the flux quantum {X g S> £o) • 

The analytical results obtained below allow a simple qualitative interpretation: after 
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cooling in a field below the glass transition temperature T g the system is locally trapped in 
one of the many metastable states of the formed glassy phase (Fig. 1). The applied current 
plays the role of the force acting on a particle in this relief, since the curvature of each valley 
is finite, the small applied force drags the particle only a little up the slope increasing its 
energy; it translates to a non-zero p s of a superconducting network. When the dragging 
force exceeds a maximal slope of the valley the particle becomes unstable and crushes to 
the next valley; for the superconductor it means that the state with a current exceeding the 
critical one becomes unstable, the phase slips and eventually a completely different state 
is formed. Before such instability occurs, the motion is reversible, i.e. the decrease of the 
current would lead to the same initial state, implying that the dissipation is due to these 
crushes only. Qualitatively this picture is similar to the vortex pinning and to the stick and 
slip model of dry friction. The role of the generalized force for this state is played by the 
variations of the vector potential A A (in the gauge VA = 0), not by magnetic field because 
the responses to it are non-local; in this respect the glassy state behaves as a superconductor, 
where the current obeys London equation j = —p s A. 

Neglecting the effects of field penetration within each granule and we assume that each 
granule is completely specified by the phase of the order parameter 0j. The interaction 
between granules is due to the Josephson couplings J if 

H = \ Yl Jijco^i ~ <f>j ~ a ij) (!) 

z a 

where = ^ // Adr is the phase difference induced by electromagnetic field A and $ — 
nhc/e; below we shall use the rescaled dimensionless vector potential a = 27r£ AA/<I> 
instead of A A. It is also convenient to introduce the effective couplings = J^e"*" , then 
the interaction energy expressed through the variables Si = e 1 ^ acquires a form reminiscent 
of the spin glasses: H = ^J2i,jJijS*Sj. The large background magnetic field H Tq ^> $ 
ensures that the couplings are random and correlations between different couplings are 
negligible, moreover the variation of the electromagnetic field A affects the correlation of 
the effective coupling on the same link: 
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Jij(t)J k i(f) = SikSflKfaftl - (Ao^(t) - A aij (t')Y}. (2) 

Here K(r) is a smooth function which determines the interaction range £o via £o = 
2 J K{r)r 2 d d r / J K(r)d r. In deriving (|2]) we assumed that the temporal variation of the 
induced phase difference is small: Aa^-(i) <C 1 and kept only leading terms in Aay(£), 
because, as we show below, in the vicinity of T c the phase slip instabilities start to occur 
long before Act^ (t) exceeds unity. 

Since the details of the microscopic dynamics are not important at the macroscopic time 
scales, we shall assume the simplest overdamped phase dynamics: 

i& = £ ImJijjte-M + tj sin fa. (3) 
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Here the damping coefficient F can be expressed through the resistivity R of the individual 
junction via T « where n is the number of granules per unit volume. In the equations 
(§) we have introduced an auxiliary local field rj(t) which is a convenient tool to probe 
the phase dynamics. As we show below, the current induced by the electromagnetic field 
variation can be conveniently expressed through the response function 

^[# (il )]) _ e 

Srj(t 2 ) 

to the field r\ and the correlator 

q(h,t 2 ) = (expli&ih) - 0(t 2 ))]) - C(h - t 2 ). 

Here we define A(ti,t 2 ) and q(ti,t 2 ) so that they describe only the long time dynamics of 
the phase, for this purpose we have subtracted the short time correlators G{t\ — t 2 ) and 
C{ti — t 2 ) decaying at scales \t\ — t 2 \ ~ 1/r. With this definition A(ix,t 2 ) and q(ti,t 2 ) 
are zero above the transition temperature where the long term memory is lost, below T g 
q(t, t) coincides with the Edwards- Anderson order parameter. The appearance of non-zero 
anomalous response function A(t 1? t 2 ) below T g is a hallmark of the glass state, the particular 
form of both A(t±, t 2 ) and q(ti, t 2 ) depends on the preparation of the state, i.e. on the path 
in (T, H) plane which lead to the state. Qualitatively it means that the realized metastable 



states are different for different cooling procedures leading to the same final temperature 
and field. Note, however, that within the MF approximation the state does not depend on 
the cooling rate (as long as it is slow compared to 1/T) but is sensitive only to the path in 
(T, H) plane. 

We shall find that in the model (0-^) the penetration depth corresponding to magnetic 
field variations is 

A, = P= r~^ } (4) 
Air^TcnTg^o 

where r = (T g — T)/T <C 1; the linear regime of the field penetration ends for field variations 
exceeding 



H e « v^ixnTg r 5/2 (5) 

where v = 0.26. At larger field variation the Bean state is formed with the small critical 
current 

27rc£ nT g 4 

j c = 7 — t 7 = 0.065 (6) 

The results (f|ffD are the quantitative results of this letter. Note that j c ~ cH c / \ g as usual. 

In order to obtain we shall derive a closed system of MF equations for the functions 
A(ti,t 2 ) and g(ti,t 2 ), discuss their solutions and finally express X g and j c through A(ti,t 2 ) 
and q(t 1 ,t 2 ). 

The equations for A(ti,t 2 ) and g(ti,t 2 ) for the infinite range Josephson network were 
obtained in Ref. and their analogs for the Sherrington-Kirkpatrick model in Ref. 
Here we sketch the generalization of the derivation of Ref. for the finite range model 
with spatial variations 0. We follow the short cut suggested by Kurchan and use the 
super symmetric representation of the spin glass dynamics; in this approach all correlators 
are deduced from the disorder-averaged generating functional Z = J exp(F{Q})VQ: 



F{Q} = 2nT g J d 3 r 



^IKoV + a 2 -a 1 ]Q(6 1 ,e 2 )| 2 + 



](r 1 +r 2 -q 1 -q 2 )Q 2 (Q 1 ,e 2 ) + ^-Q\Q 1 ,Q 2 )\dQ 1 dQ 2 + 
4 lo 

I J Q(e a , e 2 )Q(e 2 , e 3 )Q(e 3 , ei)deide 2 de 3 } (7) 

where we denoted r(ti) = Ti, a(ti) = a 1; g^tx) = g 1; etc. In the functional (|7|) we 
introduced Grassman coordinates 6, 9 which are supersymmetric partners of the time axis 
and used shorthand notations denoting the set (t, 6, 8) by 0. The functional F{Q} describes 
the long term dynamics, so we neglected in ([7]) the terms containing the time derivatives with 
coefficients ~ 1/T. The superfield Q{Qi,Q 2 ) has many bosonic and fermionic components, 
but the only components remaining large in the limit of slow field variations are bosonic 
components related to the correlators A(ti,t 2 ) and q(ti,t 2 ): 

Q(01, 6 2 ) = q(t u t 2 ) + Mi A(*2, h) + 6 2 6 2 A(t l5 t 2 ). 

The MF approximation is equivalent to a saddle point approximation for the generating 
functional Z, varying the field Q(0i, Q 2 ) we find the MF equations: 

\2q 2 (t u t 2 ) ~ r\ - rl - |( ai - a 2 ) 2 ]A(t 1 ,t 2 ) 
+|/dtA(t 1 ,t)A(t,t 2 ) = 



r |g 2 (ti,t 2 ) - r\ - rl - |( ai - a 2 ) 2 ]g(t 1 ,t 2 ) (()) 
+| J d*(g(* a , t) A(t, ti) + q(t, h)A(t 2 , t)) = 



Here the functions q(ti,t 2 ), A(ti,t 2 ) depend implicitly on the spatial coordinate r only 
through the a(r)-dependence; this simplification is due to the neglect of the gradient terms 
which is justified in MF approximation. The Eq. (fUD can be transformed to the slow 
cooling equations of infinite range model Ref. by a substitution a(t) — > 5H(t), A(ti,t 2 ) — > 



|A(ti,t 2 ) 



j(t) = jj ^ , in the MF approximation it simplifies to j(t) = 



A=0 



To find the current in a final state we use a general procedure : we replace the vector 
potential A by A + A66 and perform the variational derivative of the generating functional 

evaluated on the 

saddle point solution Q(Qi, Q 2 ) of equations (|Sf|5]). Evaluating the derivative we express the 
current through these solutions: 
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j = - nT 9^^ 1 9(*i > *2)A(t lj t 2 )(a(ti)-a(t 2 )) dt 2 (10) 

The Eqs. (|-|9|) and fllQf) are analogs of the Ginzburg-Landau equations for the glassy phase 
of superconductor. 

For a very weak variation of the external field we may replace A and q in this equation 
by the solution of Eqs. in constant field (a = 0); in this case the induced current is 

related to the vector potential by London equation: 

j = -pPSA (11) 



Pi = J§— £ J 9b(*i,*)Ao(*i,*)* = gg " t 3 



Here we used the exact solution (A, q) of the (||-||) for the cooling at constant field @. Thus, 
the superconductor in the glassy state shields the variation of the magnetic field; using the 
superfluid density (|TT|) we find \ g @. 

The solution of the equations (f§§|) changes completely when a becomes of the order 
of r, in this regime we must solve these non-linear equations from scratch and evaluate 
current using (p~0|) . Dimensional analysis shows that for the field variation that takes place 
at constant temperature the current is 

2nc^onT g 4 ^ fa 



j(a) = -°*^t^ T ■ T \r) (12) 

Here e a = a./\a\ and the dimensionless function T(a) was evaluated numerically, we show 
its plot in Fig 2; note that T(a) = a at small a < 1, reproducing (O). 

The current reaches the maximum value for a c = 0.18r, we were not able to find the 
stable numerical solution of the equations (fHI) for a > a c . As long as the solution of the 
equations (f|||) exists, the changes induced by the variations of a are reversible, i.e. if one 
changes the field a from its initial zero value at to to a(ti), then to the maximal value a(t m ) 
and then back to a(i/) = a(ti), the current (and all other properties) in the final state are 
the same as at t\. To prove it we define function $(t) so that a(#(£)) = a(t) and $(t) < t m 



(note that with this definition #(£) = t at t < t rn ). Suppose that A (ti,t2), qo(ti,t 2 ) are 
solutions of for the monotonic variation of the field (ti,t 2 < t m ), then the solution in 
the whole range ti , t 2 < i/ is given by 

A(ti,t 2 ) = 0(*i - t 2 )A (7?(t 1 ), #(t 2 ) 
q(t 1 ,t 2 )=q Q ($(t 1 ),$(t 2 )) 

which is verified by the substitution into (|^-0). Inserting the solution ( |T3D into (0) we see 
that j($(t)) = j(t) proving the reversibility of small variations of a < a c . 

Clearly, large variation of magnetic field should eventually lead to the current dissipation 
and a loss of reversibility. Thus, the solution of the Eqs. (§-0) can not exist beyond some 
critical value of a. This can be proved formally, the proof goes in three steps: (i) for a large 
constant variation of the magnetic field all memory of the initial couplings Jy(to) is lost (see 
(0)), so the order parameters depend only on the time difference: A(ti — t 2 ), q(ti — t 2 ). Now 
we find the asymptotics of these functions as t — > oo; (ii) we prove by the direct substitution 
into Eq. that solution which is identically zero at large times (A(£) > at t < tf,, but 
A(t) = at t > tb) is impossible; and finally (iii) we show that at t — ► oo Eq. @ does not 
admit decreasing solution (instead it gives A(£) oc t which make no sense). Unfortunately, 
this does not establish the maximal value of a which admits the solution of Eqs. (f|||), but 
it is likely that this value coincides with a c at which the current is maximal and numerical 
stability of the solution of is lost. 

The Eqs. (§|-0) were derived under the assumption that the variations of the phase vari- 
ables happens on the same large time scales as the variation of the external field. The 
absence of the solution of these equations can only mean that this assumption no longer 
holds and a fast phase slip occurs at t c when a(t c ) = a c . The state forming at t > t c has 
no correlations with the state before the catastrophe, this state resembles the nascent state 
obtained after a new fast cooling. Since the rate of cooling does not affect the forming 
state as long as this rate is slow compared with microscopic times scales, the nascent state 
can be described by the solutions of the slow cooling equations (§-§) in which the initial 
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period of slow cooling in a constant magnetic field to > t > was compressed to a infinitely 
narrow time interval (t c , t c + e). The solution of Eqs. (§-|9|) for such cooling is related to the 
solution for unit rate cooling by rescaling: A(ti, t 2 ) = A (ii — t c , T(h))^ ~ A(ti)S(t 2 — t c ) 
for t c < t 2 < t c + e where t c = t c - t , A(ii) = j to A (ti - £ c , £ 2 )G?t 2 and A (ti, £2), ?o(*i, £2) 
is the solution for the initial period (t c > t\ > and t c > t 2 > 0) of the field variation. 
Therefore the solution after the phase slip is given by 

A(ti, h) = 9(tx - t c )6(t 2 - t c )A (ti - t c , t 2 - Q+ 
Hh)6(t 2 - t c ) 

q(t u t 2 ) = 9(h - t c )6(t 2 - Qqoih - t e , t 2 - Q 
The continuing the field variation at t > t c results in a current growing as j(a— a c ) until the 
variation of a exceeds 2a c , then the second phase slip occurs and so on. 

For a very small variation of the external field the screening current falls off exponentially 
inside the sample: j = j(0) exp(— X g z). For larger fields j(0) becomes comparable with j c 
and the non-linear differential equation = ( 2?r J° A ) should be solved to find the 
current profile. At even higher fields AH > H c the current at the boundary exceeds j c 
and the dissipation begins, the field H c is somewhat analogous to a critical field of a type 
II superconductor, but here the fields {H cl ) for which vortices penetrate into the sample 
and the fields (H c ) that produce the pair breaking current coincide. The vector potential 
corresponding to H c equals $ a c/(27r£o) at the boundary and decays on scales X g , so H c w 
' ^° ^ n< ^ ^ ne numerical coefficient we solved numerically the equation for the current 
profile and got @. 

In this letter we did not consider two important physical effects present in real networks: 
(i) the magnetic field is shielded by individual granules leading to field inhomogeneity at 
small scales and (ii) exponentially rare transitions over the energy barriers missed by MF. 
It is not clear to us how large are these effects in real materials. 

Experimentally, the study of Josephson networks is complicated by the difficulty of sep- 
arating the response of the individual granules from the collective effects of Josephson cur- 
rents. Recently Ocio and Leylekian have surpassed this difficulty || preparing the state by 
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cooling the sample in a large field and then studying the responses to a very small varia- 
tions of field and temperature. They found that the reaction of this state to the change of 
magnetic field AH can be characterized by large \ g and a small, but non-vanishing, j c , e.g. 
they found a complete shielding of magnetic field variations. 

In conclusion we found that the effect of the variation AH of the external magnetic field 
on the Josephson network cooled in a constant magnetic field is qualitatively similar to the 
effect of the external field on the type II superconductor with an extremely strong pinning. 
Quantitatively, however, the superconductivity of the Josephson network is much weaker: 
the expressions for penetration depth, critical current, critical field contain large powers of 
reduced temperature and small numerical factors. For instance, ^ = -, % = ^P-^-t 5 / 2 
where A , j® and T c are the penetration depth, critical current and transition temperature 
of the same network in zero external magnetic field and we compare properties for the same 
reduced temperatures r. 

We are grateful to L. Leylekian and M. Ocio for numerous inspiring discussions and to 
P. Chandra for useful critique. This research was supported by the innovation partnership 
grant of NJ and by ISF through grant M6M000. 
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FIGURES 

FIG. 1. A typical energy relief of the glassy state; before (full line) and after (dashed line) the 
magnetic field was changed. Arrows show irreversible processes. 

FIG. 2. Induced scaled current T(a) as a function of the vector potential variations a = cl/t 
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